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Approximate Euclidean Shortest Paths amid Polygonal Obstacles 


R. Inkulu’^ Sanjiv Kapoor^ 


Abstract 

Given a set V of non-intersecting polygonal obstacles in defined with n vertices, we 
compute a sketch Q. oi V whose size is independent of n. We utilize fl to devise an algo¬ 
rithm to compute a (1 -I- e)-approximate Euclidean shortest path between two points given with 
the description of V. When V comprises of convex polygonal obstacles, we devise a (2 -|- e)- 
approximation algorithm to efficiently answer two-point Euclidean distance queries. 


1 Introduction 

Given a set R = {Pi, P 2 ,.. ■, Ph} of pairwise-disjoint simple polygonal obstacles in and two 
points s and t in the free space P{V) defined by the closure of sans the union of interior of all 
the simple polygons in V, the Euclidean shortest path finding problem asks to compute a shortest 
path between s and t that lies in P{V). This problem is well-known in computational geometry 
community. Mitchell |H] provides an extensive survey of research accomplished in determining 
shortest paths in polygonal and polytope domains to that date. In the following, we assume that 
n vertices together define h obstacles of V. 

The polygonal domain V is given as input, a priori. The three variants of the problem include: 
(i) both s and t are given as input with V, (ii) only s is given as input with V, and (iii) neither s 
nor t is given as input.The type (i) problem is a single-shot problem and involves no preprocessing. 
In a type (ii) problem, the preprocessing phase constructs a shortest path map with s as the source 
so that a shortest path between s and any given query point t can be found efficiently. In the 
third variation, which is known as a two-point shortest path query problem, V is preprocessed to 
construct data structures that facilitate in answering shortest path queries between any given pair 
of query points s and t. 

In solving type (i) or type (ii) problem, there are two fundamentally different approaches: the 
visibility graph method (see Ghosh [T^ for both the survey and details of various visibility algo¬ 
rithms) and the wavefront method. The visibility graph method |23LI181[T7] is based on constructing 
a graph G whose nodes are the vertices of the obstacles (together with s and t) and edges are the 
pairs of mutually visible vertices. Once the visibility graph G is available, a shortest path between 
s and f in G is found using Dijkstra’s algorithm. As the number of edges in the visibility graph 
is O(n^), this method has a quadratic time complexity barrier. In the wavefront based approach 
da ED tta [i3], a wavefront is expanded from s till it reaches t. The wavefront method typically 
constructs a shortest path map with respect to s so that for any query point t, a shortest path from 
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s to t can be found efficiently. Considering that the algorithms to find optimal shortest paths are 
complicated, simple approximation algorithms were devised in mm- 

The two-point shortest path query problem within a given simple polygon was addressed by 
Guibas and Hershberger mj. Their result preprocesses the given simple polygon in time 0{n) 
and constructs a data structure of size 0{n) and answers two-point shortest path distance queries 
in time O(lgn). The exact two-point shortest path queries in polygonal domain were explored 
by Chiang and Mitchell [S]. This result primarily devises two algorithms: one to construct data 
structures of size 0(n^} that facilitates in answering distance queries in 0(h-\-ign) time; the other 
algorithm builds data structures of size 0(n + h^) and yields O(hlgn) query time. In both of these 
algorithms, a shortest path itself is found in additional time 0(k), where k is the number of edges 
in the output path. Guo et al. [7] preprocesses T'(7^) in 0(n^ Ign) time to construct data structure 
of size O(n^) to answer two-point Euclidean distance queries for any given pair of query points in 
O(hlgn) time. 

Because of the difficulty of exact two-point queries in polygonal domains, various approximation 
algorithms were devised. Clarkson first made such an attempt in [6j. Chen [3] used the techniques 
from [6] in constructing a data structure of size 0(nlgn-|- j) in o(n^/^) + 0(j Ign) time to support 
0(6-|-e)-approximate two-point distance queries in 0{}^ + time. Arikati et al. [2] gave a family 
of results with trade-offs among the approximation factor, preprocessing time, storage space and 
the query time. Agarwal et al. [T] computes an approximate shortest path in 0{n + lg(7)) time 

when the obstacles are convex. 

All through this document, to distinguish graph vertices from the vertices of the polygonal 
domain, we refer to vertices of graph as nodes. The Euclidean distance between any two points p 
and q is denoted with ||pg||. The obstacle avoiding geodesic distance between any two points p,q 
amid a set Q of obstacles is denoted with distQ{p,q). We denote the convex hull of a set R of 
points or a simple polygon R with CH{R). For any set of polygonal obstacles Q, the free space 
resultant from the closure of sans the union of interior of all the polygons in Q is denoted with 
Let r' and r" be two rays with origin at p, and respectively make 9' and 6" counterclockwise 
angles with the positive x-axis in a coordinate system. Let and ul be the respective unit vectors 
corresponding to rays r' and r". A cone Cp{r\r") is the set of points defined by rays r' and r" 
such that a point q G Cp{R,r") if and only if q can be expressed as a convex combination of 
vectors vi and V 2 - When the rays are obvious from the context, we denote the cone with Cp. 
The counterclockwise angle from the positive x-axis to the line that bisects the cone angle of Cp is 
termed as the orientation of the cone Cp. 

Our contribution 

We compute a sketch Q comprising of simple polygonal obstacles from the set V of input simple 
polygonal obstacles. The sketch has h simple polygonal obstacles: for 1 < i < h, simple 
polygon Pi £ V is approximated with another simple polygon Qi £ Q such that Qi C CH{Pi). In 
computing sketch oiV, we identify coreset Si of vertices of each polygon Pi £ V. Further, we 
form a corepolygon Qi £ fl, which is again a simple polygon, by joining every successive pairs of 
vertices of Si that occur while traversing Pi in counterclockwise order with a line segment. Note 
that when Pj is convex, the corresponding corepolygon Qi obtained through this procedure is both 
convex as well as Qi Q, Pi- Like in [T], combinatorial complexity of Q is independent of the number 
of vertices defining P. For two points s,t £ P'(P), we compute an approximate Euclidean path 
between s and t in P(n) using an algorithm that is a variant of [6]. From this path, we compute 
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another path in J^('P) and show that this new path is a (1 + e)-approximate shortest path between 
s and t amid polygons in V. The main contributions of this paper are summarized below: 

* When V comprises of convex polygonal obstacles, we preprocess those obstacles in 0{n + 

^ Ig Ig^ ~^) to construct data structures of size 0{^=) to answer two-point 

(2 -I- e)-approximate geodesic Euclidean distance (length) queries in 0(^j^lg^;^^) time, 
where e" = - 1. 

* As part of devising the above algorithm (i.e., again when V comprises of convex polygonal 
obstacles), we devise an algorithm to compute a (1 -|- e)-approximate Euclidean distance 
between two points s, t G -T('P), given the description of V, in 0{n + Ig time, where 

e' = Vl + e — 1. Eurther, a (1 -t- e)-approximate Euclidean shortest path is computed in 
additional 0{hlgn) time. 

* Further, we extend the above algorithm for computing a (1-|-e)-approximate Euclidean short¬ 
est path between s and t in 0(n -|- h{{lgn) + (Ig/i)^"'''^ -|- Ig ^)) time when V comprises of 

simple polygons. Here e' equals to \/T+~e — 1, and <5 is a small positive constant due to time 
involved in triangulating -E('P) using [15]. 

In addition, our algorithm to compute the coreset of simple polygons and obtaining a sketch of V 
may be of independent interest in devising efficient approximation algorithms to other geometric 
optimization problems. 

As mentioned above, Agarwal et al. PP devises a (1 -|- e)-approximation algorithm to compute a 
s-t Euclidean shortest path in 0(n+ -^Ig^) time when obstacles in V are convex. Our algorithm 
is applicable even when polygonal obstacles in V are non-convex. Our approach of computing 
coresets is quite different from [T]. The key differences between our algorithm and [T] are detailed 
below. For two-point shortest paths, Chiang and Mitchell |5| outputs an optimal shortest path; 
our result answers approximate two-point distance queries with time-space trade-offs with respect 
to [2j. Moreover, our result trades-off with approximation algorithms devised in [21EI- 

In our algorithm is very different from [Tj.Let the polygonal domain V be defined with 
convex polygons Pi, P 2 , ■ ■ ■, Ph- In this algorithm as well as in [T], Pi is approximated with Qi, for 
every 1 < i < h. However, for every 1 < i < h, in our algorithm Qi C Pi whereas in P, Pi C Qi. 

Let the new polygonal domain H be defined with simple polygons Qi, Q 2 , ■ ■ ■, Qh- Unlike ra , in 
computing H, our algorithm does not require using plane-sweep algorithm to find pairwise vertically 
visible simple polygons of V. Our algorithm partitions the boundary of each convex polygon P into 
a set of contiguous patches so that the angle subtended by any two points belonging to same patch 
is upper bounded as a function of e. For any two points s,t G -T('P), our algorithm guarantees 
that distn{s,t) < dist'p{s,t) < (I -|- e)dist^{s,t). For any two points s,t G J^(H), algorithm in 
P guarantees that distQ{s,t) < (I -|- e)dist-p{s,t). To find a shortest path amid H, our algorithm 
does not use algorithm from |I2| . Instead, our algorithm achieves the said approximation using 
the spanner constructed from cone Voronoi diagrams (CVDs) [6|. Apart from computing sketch 
H of P, as compared with [H], number of cones per obstacle that participate in computing CVDs 
amid P{^) is further optimized to achieve the above mentioned time and space complexities. This 
is achieved by exploiting the convexity of obstacles together with the properties of shortest paths 
amid convex obstacles. Further, in our algorithm, for any maximal line segment with endpoints 
P, r” along the computed (approximate) shortest path p amid obstacles in H, if r'r" lies in Pi — Qi 
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for any 1 < i < /i, then before outputting p we replace line segment r'r" with geodesic shortest 
path between r' and r” in J-{V) i.e., with a shortest path between r' and r” along the boundary of 

Pi. 

More importantly, the definition of coresets and corepolygons of convex polygonal obstacles is 
extended to simple polygons using the decomposition of P{P) into hourglasses [HI [IHl [H] . And, 
the sketch n is defined to comprise simple polygons, like in the convex polygonal case each Qi € Q 
correspond to an obstacle Pi G P and is obtained from the coreset Si of Pi. The scheme designed 
in Agarwal et al. [T] does not appear to extend easily to the case of simple polygons as they use the 
critical step of computing partitioning planes between pairs of convex polygonal obstacles from V. 

The algorithm for single-shot approximate shortest path computation when obstacles are convex 
polygonal is described in Section Section describes algorithm for computing an approximate 
Euclidean shortest path amid simple polygonal obstacles. The preprocess-query algorithm to answer 
approximate Euclidean distance queries is described in Section]^ The conclusions are in Section 

2 Approximate shortest path amid convex polygons 

In this Section, we suppose that every polygon belonging to P is convex. We use the following 
notation from Yao [H]. Let C be the set of cones with disjoint interiors partitioning where each 
cone has an apex at the origin and the cone angle is upper bounded by a value that is a function 
of e. Each cone that we refer in this paper is a translated copy of some cone in C. When a cone 
C G C is translated to have apex at a point p, the translated cone is denoted with Cp. We show 
that 0(^^) vertices (a dehned later), selected from the set of vertices defining P, together with 
the select set of cones introduced at these vertices, suffice to compute an approximate shortest path 
between any two points in P(P) with the desired accuracy. 

2.1 Computing a sketch of V 

First, we detail an algorithm to compute a coreset of P. Let ej, e^+i,..., be a sequence C of 
edges when the boundary of a convex polygon Pi is traversed in counterclockwise order such that 
for every j < I < k — 1, vertex vi is common to edges ei and e^+i. The absolute value of difference in 
angle made by Cj and with the positive x-axis is defined as the angle subtended by C. Let Ilj be 
a partition of the boundary of a convex polygon Pi into a collection of [\ contiguous sections, 
patches, such that the angle subtended by any contiguous section is upper bounded by with 
a defined in terms of e later. 

Lemma 2.1 For any two points p and q that belong to any patch R Glij, the geodesic Euclidean 
distance between p and q along R is upper bounded by (1 -|- ae)||pg'||. 

Proof: Let P be the edge on which p lies and let e” be the edge on which q lies. Let c be the point 
of intersection of normal to e' at p and the normal to e” at q. Since p and q belong to the same 
patch, the angle between cp and cq is upper bounded by ^/ae. The geodesic distance between p and 
q along R is upper bounded by < (1 + f)\\PQ\\ < (1 + ae)\\pq\\. □ 

For each obstacle Pi, the coreset Si of Pi comprises of two vertices chosen from each patch in 
Ilj. In particular, for any patch tt G Ilj defined with the sequence Vj,Vj+i,... ,Vk-i,Vk of vertices 
along an obstacle Pi, both the vertices Vj and Vk belong to the coreset Si of Pi. The coreset S of 
P is then simply IJ^ Si. 
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Observation 1 The size of coreset S ofV is 0(^J=). 

Noting that the complexity of CH{Si) is upper bounded by the complexity of Pi, to achieve 
the efficiency, for every 1 < i < h, our algorithm uses corepolygon Qi = CH{Si) in place of Pi. 
Let n be a set comprising of corepolygons corresponding to each of the polygons in V. Then is 
the sketch ofV. The following Lemmas show that we can achieve a (1 + e)-approximation amid V 
while using the sketch of V. 


Lemma 2.2 Letv',v'' he any two vertices of obstacles inQ. Then, dist-p{v',v”) < {l+ae)distn{v', v" 


Proof: Let vi and V 2 be any two successive vertices along a shortest path between v' and v" amid fl. 
Let O be the set of obstacles intersected by line segment viV 2 in V. Suppose that vi and V 2 belong 
to obstacles Pj and Pk respectively. Since the line segment viV 2 does not intersect with the interior 
of the CH{Sj) or CH{Sk), it intersects with at most one patch belonging to set Ilj of patches of Pj 
and at most one patch belonging to set 11^ of patches of Pk- (Refer Fig. [^) Let the line segment 
viV 2 intersect with a patch R G lij at a points vi and r. Then from Lemma 2.1, the geodesic 


distance between vi and r along R is upper bounded by (1 + ae)||uir||. Analogously, let the line 
segment viV 2 intersect with a patch R' G 11^ at points V 2 and r'. Then the geodesic distance from 
V2 and r' is upper bounded by (1 + ae)||t'2r'||. For any convex polygonal obstacle Pi in O distinct 
from Pj and Pk, let p',p” be the points of intersection of viV 2 with the boundary of Pi. Since the 
line segment viV 2 does not intersect with the interior of the convex hull of coreset corresponding 
to Pi, both p' and p'' belong to the same patch, say R” G 11;. Then again from Lemma 2.1 
geodesic distance from p' and p" along patch i?" is upper bounded by (1 + ae)||pV'||. 


the 



Figure 1: A line segment viV 2 of a shortest path amid SI intersecting three patches belonging to obstacles in V 


Let TTi, 7r2,..., TTfc be the set IT of patches intersected by viV 2 , and let p(,pf be the points of 
intersections of viV2 with any patch tt; G LI. Then dist-p{p'-,p”) added with llpfPi+ill 

is upper bounded by (1 + Q!e)||uit; 2 ||- Let vi,... ,Vr be the vertices of P that occur in that or¬ 
der along a geodesic shortest path between vertices v',v” G P amid 11. Then distp{vi,Vr) = 
YllZl dist-p{vi,Vi+i) < (1 + ae) distQ,{vi, Uj+i). □ 


Since P{P) C 7^(11), every path amid convex polygonal obstacles in P is also a path amid 
convex polygonal obstacles in H. This observation leads to the following. 

Lemma 2.3 For any two vertices v',v" ofP, distQ{v',v") < dist-piv',v''). 

Thus we ensure that for source and destination vertices s,t G P, shortest path computed amid 
n achieves (1 -|- ae)-approximation. 

Lemma 2.4 Let P be a collection of h convex polygons in with n vertices and let s and t 
be two points in P{P). Then the sketch S of P with cardinality 0{^^) suffices to compute a 
(1 -|- ae)-approximate shortest path between s and t in P{P). 
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Proof: Immediate from Observation!^ Lemma 2.2, and Lemma 2.3 


□ 


2.2 Computing an approximate shortest path using the sketch VL o^ V 

Our algorithm relies on [6]; hence, we give a brief overview of the algorithm from [6] to construct a 
Euclidean spanner Giy, E) from the given set of obstacles V. Noting that the endpoints of segments 
of a shortest path in comprises of only vertices of V, the node set V is defined as the vertex set 
of P. Let C be the set of 0(j) cones apexed at origin that together partition E{V). Let C £ C 
be a cone with orientation 9 and let C £ C he the cone with orientation —9. A cone Voronoi 
diagram CVD is constructed corresponding to each cone C £ C such that for a cone C £ C and set 
K of points, we let CVD[C,K) be a partition of PiV) where for every point p £ K there is an 
associated region Rp C py)^Rp £ CVD{C, K). Rp is defined using the cones in C. Indeed, the 
region Rp comprises of all points q £ P{P) such that p is the closest vertex in C' among points 
in K. For a given cone, Cy, let V' be the set of vertices of V that are visible from v and that lie 
within the cone C^. The vertex in V that is closest to v, say v', is said to be the closest vertex in 

to V. For every vertex v oiV and for every cone Cy, an edge e joining v and a closest vertex in 
Cv to V, say v', is introduced in E with its weight equal to the Euclidean distance between v and 
v'. Apart from these, no additional edges are added to E. The result in [6] proves that if d is the 
obstacle avoiding geodesic Euclidean distance between any two vertices, say v' and v", of P, then 
the distance between the corresponding nodes P and v” in G is upper bounded by (l + ejd. Further, 
|6] computes GVD{G, K) using plane-sweep in 0(| ATj Ig |iL|) time; and, well-known planar point 
location structures facilitate locating any point q in any GVD while help in designing shortest path 
query algorithms. 

By limiting the number of vertices of P at which the cones are initiated to coreset S of vertices, 
our algorithm improves the space complexity of the algorithm in [6]. Further, by exploiting the 
convexity of obstacles, we introduce 0{^^) cones per obstacle, each with cone angle O(y^), and 
show that these are sufficient to achieve the claimed approximation factor. 

Let uo, ui,..., Ufc be vertices such that ui,..., Vk-i belong to an obstacle P and vo,Vk belong 
to obstacles P' and P" respectively, for P p P' p P". Also, let vqVi .. .Vk he a subpath vr of a 
shortest path. Since P is a convex polygon, ZvjVj+iVj +2 interior to P is less than tt radians for 
every j £ [0, A; — 2], i.e. the subpath is convex w.r.t. P at each of the vertices^vi,V 2 , ■ ■ ■ Vk-i- 

Let u be a vertex of P that belongs to core¬ 
set Si of convex polygon Pi. Let P,v, v" be the 
vertices that respectively occur while traversing 
the boundary of Pi in counterclockwise order. 

Also, let C' be the cone defined by the pair of 
rays (uu', —vv^') and let G” be the cone defined 
by the pair of rays {vP, —vP). For a coreset 
vertex v £ S, a cone G £ C is said to be ad¬ 
missible at V whenever n C" or Gy n G" is 

non-empty. (See Fig. i) Note that any shortest path is convex at v with respect to Pi. Thus if q 
is any point in P{P) such that q is not visible to p amid P and a shortest path between p and q 
passes through vertex v of Pi, then there exists a shortest path from p to q such that one of its line 
segment lies in C' and another line segment of that path lies in G". Hence, in computing geodesic 
shortest path amid P, it suffices to consider admissible cones at vertices of P. 

Note that whenever two points s and t between which we intend to find a shortest path are 


Cv 



C" 


C 


Figure 2: Illustrating an admissible cone Cv incident to 
a coreset vertex v of an obstacle 
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visible to each other, the line segment st needs to be computed. To facilitate this, for every 
degenerate point obstacle p, every cone C with apex p is considered to be an admissible cone. 

The same properties carry over to the domain 0 as well. For any two points pi and p 2 in T'(O), 
suppose Pi and p 2 are not visible to each other. Consider any shortest path tt between pi and p 2 - 
For any line segment ab in vr, ab is either an edge of O or it is a tangent to an obstacle O. In the 
latter case, ab belongs to an admissible cone of O. 

The following Lemma upper bounds the number of cones introduced. 


Lemma 2.5 The number of cones introduced at all the obstacles of Tl together is 

Proof: Let be a ray with origin of the coordinate axes as its endpoint. (See Fig. [^) For any 
two distinct vertices v' and v" of a convex polygon P, let be the ray parallel to with origin 
at v' and pointing in the same direction as and let r\yi be the ray parallel to r with origin at v” 
and point in the same direction as r. Also, let v[ precede v' (resp. v'{ precede v”) and succeed 
v' (resp. v'f succeed v") while traversing P in counterclockwise order. Since P is a convex polygon, 
if FJ/ belong to cone defined by v\ v' and then rf^/ is guaranteed to not belong to cone defined 
by v''v 2 and 



Figure 3: Illustrating that a ray parallel to r can exist in only one admissible cone per obstacle 
Extending this argument, if a cone Cy' is admissible at v' then the cone cannot be ad¬ 
missible at v". Since the number of coreset vertices per obstacle is the number of cones 

introduced per obstacle is 0{^^). Further, since there are h convex polygonal obstacles, number 
of cones at all the obstacle vertices together is 0{-i=). □ 



Figure 4: Illustrating an edge of the spanner 

Next, we describe computing the spanner G{V = 5 U S',E). The set S comprises of nodes 
corresponding to coreset S. For every v £ S and for every admissible cone Cy, let V' be the set of 
points on the boundaries of obstacles of Q that are visible from v and lie in Cy. (See Fig. 0) The 
point in V' that is closest to v, say p, the closest Steiner point in Cy to v is found and p is added 
to Sh An edge e between v and p is introduced in E while the Euclidean distance from v and v' is 
set as the weight of e in C. Let p be located on a convex polygonal obstacle P. Eurther, for every 
Steiner point p, let v' (resp. v") be the coreset vertex or Steiner point that he on the boundary of 
P and occurs before (resp. after) p while traversing the boundary of P in counterclockwise order. 
Then an edge P (resp. e'') between p and v' (resp. p and v"') is introduced in E while the Euclidean 
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distance from p to v' (resp. p and v") is set as the weight of e' (resp. e") in G. The set S' comprises 
of all such Steiner points. Note that both \V\ and \E\ are 0(-P=). 


Lemma 2.6 Let G be the spanner constructed from Ll. Let distcip',p") be the distance between p' 
and p" in G. Then for any two points s,t G distQ{s,t) < distG{s,t) < (1 + a€)distQ{s,t). 


-1 


Proof: Theorem 2.5 of [6] concludes that to achieve (1 + ae)-approximation, sin'0 — cosip < 

Expanding sine and cosine functions for the first few terms yield —1 + 1 ^ + ^ < Solving 

the quadratic equation in ip yields ip < ^fae. Hence, the choice of the cone angle of cones and the 


cardinality of C in Lemma 2.5 


We claim that introducing a subset of cones (admissible cones) rather than all the cones as used 
in does not affect the correctness. Let p and q be vertices of convex polygons Gp and Gg respec¬ 
tively. Suppose pq is a line segment belonging to a shortest path in the spanner computed in [6]. If 
pq is a tangent to convex polygonal obstacle Pi at p, then pq belong to an admissible cone Gp at p. 
Similarly, if pq is a tangent to convex polygonal obstacle Pj at q, then pq belong to an admissible 
cone at q. Otherwise, there exists a line segment in admissible cone apexed either at a vertex of 
Pi or at a vertex of Pj which yield a shorter path from source s to g without using pq line segment. □ 

Once we find a shortest path SPq, amid obstacles in LI using spanner G, as for the proof of 


Lemma 2.3, we transform it to a shortest path amid obstacles in V. Since there are 0{h) obstacles in 


Ll, SPq contains 0{h) tangents. Let this set of tangents be T. We need to find points of intersection 
of convex polygons in V with line segments in T. Whenever a line segment I G T and a convex 
polygon Pi G V intersect, say at points p' and p", we replace the line segment between p' and p" 
with the geodesic shortest path between p' and p" along the boundary of Pi. Analogously, for every 
line segment I G SPq — P belonging to an obstacle Pj G LI, we replace I with the corresponding 
geodesic path along the boundary of Pj. We use the plane-sweep technique [8] to find the points 
of intersections of line segments in T with the convex obstacles in V. 

As part of plane-sweep, a vertical line is swept from left-to-right in the plane. Let L (resp. R) 
be the set of leftmost (resp. rightmost) vertices of convex polygons in V. Initially, points in L and 
R together with the two endpoints of every line segment in P are inserted into the priority queue, 
say Q. The event points are scheduled from Q using their distance from where the sweep line is 
initially placed as priority. As the algorithm progress, the event points corresponding to L, R, and 
the endpoints of line segments in P are deleted from Q. The algorithm terminates whenever the 
Q is empty. As described below, the intersection points between the line segments in P and the 
convex polygons in P are added to Q with the traversal of sweep line. The sweep line status is 
maintained as a balanced binary search tree B. We insert (resp. delete) a line segment in P or 
convex polygon in P, say s, to B whenever leftmost (resp. rightmost) endpoint of s is popped 
from Q. Since before an intersection occurs between a line segment I from P and a convex polygon 
P from P, it is guaranteed that I and P occur adjacent along the sweep line, we update event- 
point schedule with an intersection between I and P whenever I and P are adjacent in the sweep 
line status. By using the algorithm from Dobkin et al. [9], we compute the possible intersection 
between I and P. If they do intersect, we push the leftmost point of their intersection to Q with the 
distance from initial sweep line as the priority of that event point. Further, we store the rightmost 
intersection point between I and P with the leftmost point of intersection as satellite data. If the 
leftmost intersection point between I and P pops from Q, we compute the geodesic shortest path 
along the boundary of P between the leftmost intersection point and the corresponding rightmost 






intersection point. Further, whenever I and P become non-adjacent along the sweep line, we delete 
their leftmost point of intersection from Q. 

Theorem 2.1 Given V and two points s and t in P{P), computing a (1 + e)-approximate distance 
between s and t takes 0{n + ;^lg^) time for e' = y/1 + e — 1. Further, within an additional 
0{hlgn) time, o (1 + e)-approximate shortest path is computed. 


Proof: From Lemma [2^ we know that dist'p{s,t) < distQ{s,t) < {l + ae)dist-p{s,t). Let G be the 
spanner constructed. From Lemma 2.6, we know that distQ{s,t) < distG{s,t) < {1 +ae)distQ{s,t). 
Combining these two inequalities yields distp{s,t) < distcisp) < (1 + ae)‘^dist-p{s,t). To achieve 
(1 + e)-approximation, we set a = ^(\/l + e — 1). From here on, we denote ae with e'. 

Finding coreset vertices and computing corepolygons takes 0(n) time. The number of coreset 
vertices is O(^). The number of cones per obstacle is 0{^). Therefore, the total number of 

cones is 0{^=). For any cone C ^ C and for any corepolygon O G fl, at most a constant number 
of vertices of O are apices to cones that have the orientation of C. Considering a sweep-line in the 
orientation of C, the sweep-line algorithm to find the closest Steiner point to apex of each cone C 
(whenever an obstacle intersects with C) takes 0{h\gh) time. Hence, computing closest Steiner 
points corresponding to all the cone orientations in C together take \gh). 

The number of nodes in the spanner G is O(^). This include coreset vertices and at most 
one closest Steiner point per cone. As each cone introduces at most one edge into G, the number 
of edges in G is 0(-^). Finding a shortest path between s and f in G takes 0(-^lg-^) time. 


Hence, computing the (1 -|- e')-approximate distance between s and t takes 0{n + Ig -^) time. 

For the plane sweep, leftmost and rightmost extreme vertices of convex polygons in V are 
found in 0{n) time. There are 0(h) line segments in T, cardinality of H is 0{h), and 0(h) line 
segment-obstacle pairs (respectively from P and V) that intersect. The number of event points 
due to L,R, and endpoints of line segments in T is 0(h). If I and P become non-adjacent along 
the sweep-line, deleting their point of intersection from Q is charged to the event that caused them 
non-adjacent. The sweep-line status structure is updated at the points in R as well as with the 
rightmost endpoints of line segments in T ; there are 0(h) such event points. Analogous to the 
analysis provided for line segment intersection [8], our plane sweep algorithm takes 0(n + high) 
time. 

Due to Dobkin et al. [9], determining whether a line segment I in SPq intersects with an ob¬ 
stacle P takes O(lgn) time. The preprocessing structures corresponding to [9] take 0(n) space and 
they are constructed in 0(n) time. Further, replacing every line segment between points of inter¬ 
section with the corresponding geodesic shortest path along the boundaries of obstacles together 
takes 0(n) time altogether. □ 


3 Computing approximate shortest path in polygonal domain 

In this section, we extend the approximation method to the case of non-convex polygons. We reduce 
the problem of computing approximate Euclidean shortest path amid simple polygonal obstacles 
to that of computing the approximate Euclidean shortest path amid convex polygonal obstacles. 
This is accomplished by first decomposing P(V) into a set of corridors, funnels, hourglasses, and 
junctions HanSKle]. First, we describe these geometric structures and then detail how these help 
in the reduction. 
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For convenience, obstacles of V are assumed to be contained in a rectangle TZ. Let Tn{J^) denote 
a triangulation of -F'('P). (Refer Fig. [5l) The line segments of Tri{T) that are not obstacle edges 
are referred to as diagonals. Let G{iF) denote the dual graph of Tn{T), i.e., each node of G{T) 
corresponds to a triangle of Tri{F) and each edge connects two nodes corresponding to two triangles 
sharing a diagonal of Tri{F). Based on G{F), we compute a planar 3-regular graph, denoted by 
G^ (the degree of every node in G^ is three), possibly with loops and multi-edges, as follows. First, 
we remove each degree-one node from G{F) along with its incident edge; repeat this process until 
no degree-one node remains in the graph. Second, remove every degree-two node from G{F) and 
replace its two incident edges by a single edge; repeat this process until no degree-two node remains. 
The resulting graph is G^ (refer Fig.j^, which has 0{h) faces, nodes, and edges. Every node of G^ 
corresponds to a triangle in Tri{F), called a junction triangle. (Refer Fig. H) The removal of the 
nodes for all junction triangles from G^ results in 0{h) corridors, each of which corresponds to an 
edge of G^. This procedure leaves points s and t between which geodesic Euclidean shortest path 
needs to be computed are to be in their own corridors. 

The boundary of each corridor G consists of four parts (see Fig. [^: (1) A boundary portion 
of an obstacle Pi G V, from a point a to a point 6; (2) a diagonal of a junction triangle from h to 
a point e on an obstacle Pj G V {Pi = Pj is possible); (3) a boundary portion of the obstacle Pj 
from e to a point /; (4) a diagonal of a junction triangle from / to a. The corridor C is a simple 
polygon. Let 7r(a, 6) (resp., 7r{e, f)) be the Euclidean shortest path from a to 6 (resp., e to /) in 
G. The region Hq bounded by 7r(a, 6), 7r(e,/), be, and fa is called an hourglass, which is open if 
7r(a, b) n 7r(e, /) = 0 and closed otherwise. (Refer Fig. [^) If He is open, then both 7r(a, b) and 
7r(e,/) are convex polygonal chains and are called the sides of He', otherwise. He consists of two 
funnels and a path vre = 7r(a, b) n 7r(e, /) joining the two apices of the two funnels, and vre is called 
the corridor path of G. The paths 7r{b, x),7r{e,x),TT{a,y), and vr(/, y) are termed sides of funnels 
of hourglass He- The sides of funnels are convex polygonal chains. For any obstacle Pj £ V, let 
7r{b,x) be a side S' of a funnel F of He such that 6 is a vertex of Pj. Let 7r{b,x') be a maximal 
subpath of 7r{b,x) such that x' is a vertex of Pj. Then x' is termed a pseudo-apex of side S of 
funnel F. Note that pseudo-apex of a side S could be same as the apex of the funnel F. Let x', x" 
be two pseudo-apices of a closed hourglass He such that x' and x” are vertices of Pj. The shortest 
path between x' and x” along the boundary of He is the corridor path between pseudo-apices. 



Figure 5: Illustrating a triangulation of the free space 
among two obstacles and the corridors (indicated by red 
solid curves). There are two junction triangles marked 
by a large dot inside each of them, connected by three 
solid (red) curves. Removing the two junction triangles 
results in three corridors. 




Figure 6: Illustrating an open hourglass (left) and a closed 
hourglass (right) with a corridor path connecting the apices x 
and y of the two funnels. The dashed segments are diagonals. 


For any obstacle Pj £ V, let Cj be the union of the following: set of open hourglass sides whose 
endpoints are incident to Pj , set of maximal sections of sides of funnels whose endpoints are incident 
to Pj, and the set of line segments each of which joins pseudo-apices that correspond to He. The 
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closed and bounded polygon defined by polygonal chains in Cj contains Pj] hence, we term this new 
polygon as the expanded polygon of Pj. If every chain in Cj is convex, then IJceC- ^ CH{Pj). 
Since s and t lie exterior to closed hourglasses, as described below, in forming the spanner graph 
G, we handle corridor paths between pseudo-apices as a special case and introduce corresponding 
edges into G. 

For every convex polygonal chain C of every open hourglass that correspond to any polygonal 
obstacle Pj, endpoints of G are chosen to be in the coreset Sj corresponding to Pj. Further, similar 
to the case of convex polygonal obstacles, we compute a coreset of every such convex polygonal 
chain G and replace G with the convex polygonal chain G' that passes through those coreset of 
vertices. Essentially, we partition G into patches with angle subtended by each patch to be same as 
in the case of convex polygonal obstacles. (For details, refer Section i) For every closed hourglass 
He that has a vertex common to Pj , both the pseudo-apices of He that incident to Pj are included 
into coreset of Pj. Similar to convex polygonal chains of open hourglasses, coresets corresponding 
to maximal sections of sides of funnels whose endpoints incident to Pj are computed; and, convex 
polygonal chain that passes through these coreset of vertices of G is used. Excluding line segments 
joining successive pseudo-apices, the convex polygonal chains that pass through coreset vertices of 
each polygonal chain of Pj bound a closed region termed corepolygon Pj of Pj. Since there are 0{h) 
open and closed corridors together and since there are 0(1) convex sides per corridor the following 
is immediate. 


Observation 2 


The size of coreset S of V is 


Similar to convex polygonal obstacle case, we construct a spanner G{V, E) that correspond 
to corepolygons of V using CVDs. It is immediate to note that Clarkson’s method extends to 
corepolygons defined as above. And, for each line segment that joins two pseudo-apices along a 
corepolygon, an edge e{v',v'') is introduced into G: here, v' and v" are the nodes corresponding to 
endpoints of G; and, the weight of e is set as the distance between v' and v” along G, which is the 
length of G. Since both the pseudo-apices are included in the coreset of the corresponding obstacle, 
v',v" G V. Eor a shortest path p between any two nodes of G, before outputting p, for every edge 
e G p if both the endpoints of e correspond to pseudo-apices a', a” of a closed hourglass then we 
replace p with the shortest path along the corridor path between a' and a". Computing hourglasses 
of P{V) using [HIITHKTG] takes 0{n + h{\ghY^^ + h\gn) time (where d is a small positive constant 


due to time involved in triangulating PifP) using [E]). Extending the proof of Theorem 2.1 leads 
to the following. 


Theorem 3.1 Given a set V of polygonal obstacles and two points s and t in PifP), computing a 
{I+e)-approximate Euclidean shortest path between s andt in 0{n+h{f\gn) + {\ghY^^+{^\g^))) 

time. Here e' is Vl + e — 1, and 6 is a small positive constant. 


4 Two-point approximate distance queries amid convex polygons 

Our query algorithm constructs an auxiliary graph from the spanner network computed during the 
preprocessing phase of the algorithm. Like in the previous section, our preprocessing algorithm 
relies on [6]. We compute the approximate distance between the two query points using a shortest 
path finding algorithm in the auxiliary graph. 
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4.1 Preprocessing 


The graph G constructed as part of preprocessing in Section |2.2| is useful in finding a geodesic 
shortest path between any two vertices in V. However, since the time involved in finding a shortest 
path between nodes in G, or the space to save all-pair-shortest-paths in G are resource intensive, we 
compute a planar graph G^\V, from GiV, E) using the result from Chew |Tj. Chew’s algorithm 
finds a set E^^ C E in 0(|C| Ig |H|) time so that the distance between any two nodes of Gpi is a 
2-approximation of the distance between the corresponding nodes in G. We use the algorithm from 
Kawarabayashi et al. [19] to efficiently answer (1 + e)-approximate distance (length) queries in G^K 
More specifically, [l9| takes 0(|H| Ig^ |H|) time to construct a data structure of size 0(|H|) so that 
(2 + e)-approximate shortest distance queries are answered in time. 

Lemma 4.1 Let G be the spanner computed for the polygonal domain Ll using the algorithm men¬ 
tioned in Subsection 2.2. Let s, t be two points in E{V). Let be the planar graph constructed from 
G using W- Further, let distK{s,t) be the distance between s and t in G^^ computed using the algo¬ 
rithm from m With appropriately chosen parameters, dist'p{s,t) < distK{s,t) < {2+e)dist'p{s,t). 


Proof: From Lemma 2.4, we know that dist-p{s,t) < dist^{s,t) < (1 + ae)dist-p{s,t). Let 
distc{s,t) be the distance in G between nodes s and t of G. From Lemma |2.6l we know that 


distQ{s,t) < distG{s,t) < {l-\-ae)distQ{s,t). Let distQpi{s,t) be the distance in between nodes 
s and t of G'^K From |1|, distG{s,t) < distGpi{s,t) < 2distG{s,t). Then, as mentioned above, 
distQpi{s,t) < distK{s,t) < {1 -\- ae)distGPi{s,t). 

Combining these, dfsf'p(s, f) < distK{s,t) < {l-\-ae)^{2)dist'p{s,t). To achieve (2-|-e)-approximation, 
we set a to “ !)• Since a < — 1 for e < 1, algorithm yields a (2 -|- e)-approximation. □ 


In the following, we denote ae = — 1 with e". From here on, we suppose that there are 

0{^^) cones in C, each cone with a cone angle 0{^/G). It remains to describe data structures that 
need to be constructed during the preprocessing phase so as to obtain the closest vertex of query 
point s (resp. t) in a given cone Cg (resp. Gt). The significance of the same is explained later. To 
efficiently determine all these 0{^^) neighbors to s and t during query time, we construct a set 
of 0{^^) cone Voronoi diagrams {GVDs), each of which correspond to a cone in C. For any cone 
G € C with orientation fj, the CVD in that orientation is denoted with GVD.,p. For each G £ C, 
we consider a sweep-line orthogonal to the orientation of G. The sweep-line algorithm details are 
same as mentioned in Subsection 12.21 

Lemma 4.2 The preprocessing phase takes 0{n -\- -^Ig -\- Ig^ ■^=) time. The space com¬ 
plexity of the data structures constructed by the end of preprocessing phase is 0{^^). Here, e" is 

(2±^)V3_i. 

Proof: Computing sketch H of P takes 0(n -|- ;^ 7 ^) time. The number of cones in all GVDs is 
0{^=). Constructing spanner G involve computing CVDs, together it takes 0{^=p^= Ig time 
to compute G. Due to |1|, computing planar graph G^^ with nodes takes 

time. Computing space-efficient data structures using nni take 0(^^lg^^) time. Hence, the 
preprocessing phase takes 0{n + ^ Ig Ig^ time. 

Further, data structures constructed using m by the end of preprocessing phase occupy 0{^^) 
space. Using Kirkpatrick’s point location [20|, data structures for planar point location are of 
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0 {^=) space complexity. 


□ 


4.2 Shortest distance query processing 


The query algorithm finds the geodesic distance between any two given points s,t € J'i'P)- We 
construct a graph Gst from by introducing query points s and t as nodes in Ggt- For every 
G € C, whenever Gs (the cone with apex s) intersects an obstacle, we determine a point p on an 
obstacle, say Pi, that is closest to s in Cg. We add an edge between s and p in Gst with weight 
equal to the Euclidean distance between s and p. Let p lie between two successive coreset vertices, 
say Pi and p 2 of Pi. We add edges between p and pi and between p and p 2 with the respective 
geodesic distances along the boundary of obstacle Pi. The set Vs comprises of all such neighbors 
of s (like p) whereas each such p correspond to a cone G € C. Analogously, we define neighbors to 
t in Gst- 

The node set of Gst is 14U 14U{s, t}. The edges of this graph are of three kinds; {s} x Vs, VsXVt 
and {t} X Vt. For every edge (s, s') with s' G Vs, the weight of edge (s, s') is as defined above. The 
same is true with the edges of type {t, t') with t' € Vt. Note that every node u G 14 U 14 correspond 
to some vertex in coreset S. For every edge {s',t') with s' G Vg and t' G 14, the weight of {s',t') is 
the (2 + e)-approximate distance between s' and t'. These weights are obtained from data structures 
maintained due to |19j . 

We apply Fredman-Tarjan algorithm to find a shortest path between s and t in Gst- It follows 
from the above discussion, that this distance is a (2 + e)-approximate distance from s to t amid V. 
Computing shortest path between s and t in this way effectively equivalent to having both s and t 


available before preprocessing V. Therefore, the following is immediate from Lemma 4.1 


Lemma 4.3 The shortest path obtained between s and t in Gst is a {2 + e)-approximate shortest 
path between s and t amid set V of convex polygons. 


Lemma 4.4 The time to query the (2 + e)-approximate distance between two vertices in V is 


Proof: Locating a query point s (resp. t) in any one CVD,p takes 0(lg;^^) time. Since there 
are 0{^^) CVDs and since finding a closest Steiner point (or, vertex) in each such CVD takes 
0(lg;^) time, computing all the closest Steiner points (or, vertices) of s and t together take 
Ig '^^)- The number of nodes and the number of edges of Gst are respectively 0{^^) and 
O(^). Using [T^, finding the distance between any two nodes in take 0(^^^lg^ ^J=) time. 

There are 0 (^ 7 ) pairs of neighbors of s and t that needs to be considered in computing approxi¬ 
mate distances in G^K Fredman-Tarjan algorithm to compute a shortest path in Gst from s to t 
takes 0{^ Ig* :^^)- The Theorem follows when these asymptotic time complexities are added up. □ 


Theorem 4.1 Given a set V of h pairwise-disjoint convex polygons of total complexity n in 
preprocess V in 0{n + ^ Ig Ig^ '^fj) to build a data structure of size For 

any two query points s,t ^ -F(’F’), a (2 -|- e)-approximate distance is found in 0(^p^ Ig^ time. 
Here, e" = (I^)^/^ - 1. 
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5 Conclusions 


We have presented an algorithm to compute a (1 + e)-approximate Euclidean shortest path amid 
polygonal obstacles in The time and space complexities of our algorithm better the algorithms 
devised in Chen [3]. We also devised a (2 + e)-approximation algorithm to answer two-point 
Euclidean distance queries. Our result trades-off with approximation algorithms devised in [21E]. 
These algorithms rely on computing a sketch of the polygonal domain by extracting coresets of 
polygonal obstacles and further computing corepolygons from each of these coresets. 
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